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f^i [ Abstract. The dynamical parameters conventionally used to specify the orbit of 

1} ' a test particle in Kerr spacetime are the energy E, the axial component of the 

angular momentum, L z , and Carter's constant Q. These parameters are obtained by 



PL, 



solving the Hamilton-Jacobi equation for the dynamical problem of geodesic motion. 
Employing the action-angle variable formalism, on the other hand, yields a different 
set of constants of motion, namely, the fundamental frequencies uj r , ujg and lj^ 
associated with the radial, polar and azimuthal components of orbital motion. These 
■ frequencies, naturally, determine the time scales of orbital motion and, furthermore, 

\ the instantaneous gravitational wave spectrum in the adiabatic approximation. In 

04 ' this article, it is shown that the fundamental frequencies are geometric invariants and 

explicit formulas in terms of quadratures are derived. The numerical evaluation of 
| these formulas in the case of a rapidly rotating black hole illustrates the behaviour 

""^ ■ of the fundamental frequencies as orbital parameters such as the semi-latus rectum p, 

the eccentricity e or the inclination parameter 9- are varied. The limiting cases of 
^ ' circular, equatorial and Keplerian motion are investigated as well and it is shown that 

bJQ. known results are recovered from the general formulas. 
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1. Introduction 

In Newtonian theory, there are two basic methods of analysing the motion of a test 
mass under the attraction of a central, spherically symmetric gravitational field: Either 
Hamiltonian mechanics is applied to solve the problem in a spherical polar coordinate 
system or, for bound orbits, one may invoke a canonical transformation to action-angle 
variables. In the latter case, solving the Hamilton-Jacobi differential equation yields the 
integrals of motion and, in particular, Kepler's third law which for elliptical orbits, in 
geometric units, reads 0, [| 

/l-e 2 \ 3/2 

Mn K =(—j , (i) 
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where f^K is the angular frequency of orbital revolution, e is the eccentricity of the orbit, 
p/(l — e 2 ) its semi-major axis and M is the mass which acts as the source of the central 



A generalisation of the Hamilton- Jacobi method to the relativistic motion of a 
test mass in Kerr spacetime, as proposed by Carter [[|, shows that the motion is still 
separable and there is a complete set of constants of motion: The energy E as measured 
by an observer at spatial infinity, the axial angular momentum component L z and the 
Carter constant Q. Whereas E and L z reflect isometries of the spacetime geometry, Q 
arises from the separation of the radial and polar components of motion. The resulting 
first-order equations of motion can be used to compute the particle's worldline and, in 
limiting cases such as circular and equatorial orbits, to calculate periods of motion. 

However, no attempt has been made yet to formulate a relativistic generalisation 
of action-angle variables for geodesic motion in Kerr spacetime and to calculate the 
dynamical frequencies of arbitrary bound non-plunging orbits. Exactly these frequencies 
would be relevant to the computation of gravitational radiation emitted by stellar-mass 
compact objects in the vicinity of super-massive black holes if a perturbation theory 
based on the adiabatic approximation is employed. In this approximation one considers 
the motion as being nearly geodesic under the assumption that the time scale of radiation 
reaction, Trr, is much larger than the dynamical time scales T r , T$ and of the three 
components of orbital motion. Recently, Hughes || has raised the question what time 
scales should, in general, be chosen for T r , T 9 and T^, as there is no obvious way of 
calculating orbital periods from the integrals of motion if the orbit is neither circular 
nor equatorial. 

In the following, it is shown that suitable action-angle variables can be formulated 
within a fully relativistic framework for general orbits in Kerr spacetime and formulas 
for dynamical frequencies in terms of numerically solvable integrals are derived. 

2. Dynamical equations and integrals of motion 

Let us first review the procedure applied by Carter || in order to determine the constants 
of motion. We define the relativistic Hamiltonian for geodesic motion of a particle that 
possesses no spin in a spacetime with given metric in the absence of an electromagnetic 



In this equation, the metric components g**" are considered to be functions of the 
coordinates x a and the quantities pp are the conjugate momenta of the particle 
associated with these coordinates. Substituting the contravariant components g^ v of 
the Kerr metric in the Boyer-Lindquist coordinate representation [[?]], the explicit form 
of the Hamiltonian is given by 



field. 



field as § 




(2) 



H^ BL \x a ,pp) 



(r 2 + a 2 ) 2 - Aa 2 
2AS 



sin 2 # 



(Pt) 2 



laMr 



AS 



PtP<t> 



Celestial mechanics in Kerr spacetime 3 
A — a 2 sin 2 9. , 9 A . 1 . , 9 

where A = r 2 — 2Mr + a 2 and £ = r 2 + a 2 cos 2 0. 

A complete set of constants of motion can be determined if a canonical 
transformation $ : (x ,^) i— ► (I a , P^) is found such that the Hamiltonian becomes 
cyclic in all of the new generalised coordinates X a , i.e., (H^ BL ^ o $ _1 )(X Q ,P / 3) = 
H^ cycl \Pp), and the transformed momenta Pp are thus conserved along the worldline of 
the particle. The generator of such a canonical transformation, W(x a , Pp), is called the 
characteristic function and is determined by the Hamilton- Jacobi differential equation 
uu dWdW 2 n 

^d^ + ^ = °- (4) 
Solving this equation, yields three constants of integration, E, L z and Q, apart from 

— /i 2 /2, the value of the Hamiltonian evaluated along the worldline of the particle. 

The transformed conjugate momenta Pp can be chosen as any functions of the these 

constants, i.e., Pp = fp(—fi 2 /2,E,L z ,Q), assuming that / is bijective and C°°, and 

the generalised coordinates X a = dW/dP a are obtained from Hamilton's equations of 

motion, 

dX a OH (cycl) 

= W. = < 5) 

where each v a is a constant. The generalised coordinates are therefore given by 
X a (r) = X a (0) + v a T if the Hamiltonian is cyclic in all of these coordinates. 

For geodesic motion in Kerr spacetime, the constants of motion E := —p t and 
L z := p<j) are related to the isometries of the metric in the coordinates t and </>. E and 
L z can be interpreted, respectively, as the particle's energy and axial angular momentum 
component as seen by an observer at spatial infinity. The third constant of motion, Q, 
is obtained by separating radial and polar motion and is called Carter's constant. In 
terms of fi, E, L z and Q, the characteristic function W is given by 



where 



W = -Et + I ^ dr + / v O (W - L ,<>. ( G ) 



R = [{r 2 + a 2 )E - aL z ] 2 - A[fi 2 r 2 + (L z - aE) 2 + Q], (7) 

L 2 



e = Q 



(fi 2 - E 2 )a 2 + 



cos 2 9. (8) 



sin 2 

Since pp = gp a dW / dx a , the following set of first-order equations of motion is 
obtained |, §: 

P - a (aE sin 2 9 - L z ) , (9) 

(10) 
(11) 

L. 
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where P = E(r 2 + a 2 ) — aL z . 

If the conjugate momenta are chosen such that the Hamiltonian becomes identically 
equal to Po, i.e., H( ld \Pp) = P , then we obtain z/° = 1 and v k = 0. Substituting the 
characteristic function into the identities dW/dP = r + X°(0) and dW/dPk = X h (0) 
and adjusting the constants X a (0) in a suitable way, we find the following well-known 
integrals of motion 0, [J : 

r r' 2 J , [ 9 a 2 cos 2 6' Q , 
t-t = —dr'+ =— d9', (13) 



R 



1 f r 1 9i2 , . I T 1 dO 
*"*o =o/ ws^d^+o/ "^7^', (14) 



2j ro AVRdE 2j eo y/QdE 
1 r 1 ^'-I/^f^', (15) 



2j ro AVRdL z 2j dQ ^/QdL 
and, finally, 

In general, there are distinct turning points of the radial and polar motions which are 
asynchronously passed by the particle because r(r) and 6{t) are not periodic functions 
of time. In consequence, equation ( |i~4"D evaluates to different intervals of coordinate time 
when integrating over various cycles of radial or polar motion and, for this reason, there 
is no obvious way of calculating the time scales T r , Tg and Ta from the above integrals. 



3. The fundamental frequencies 

Even though it might seem that there is no useful notion of orbital frequencies if the 
coordinate functions representing the motion are not periodic, it is still possible to 
find a representation in which complete periodicity becomes manifest. However, such 
a representation is not found by means of a coordinate transformation in the usual 
sense, but by virtue of the more general concept of a canonical transformation of 
both spacetime coordinates and conjugate momenta. By its definition, a canonical 
transformation acts on the cotangential bundle T*(M. KcTT ) of Kerr spacetime, as opposed 
to a pure coordinate transformation between charts on the basis manifold -Mkcit- In 
analogy to phase space trajectories in non-relativistic Hamiltonian mechanics, let us 
consider the locus of all points in T*(M. Kevv ) which are passed by a particle in the 
course of its motion along a particular orbit. We shall denote the image of that set 
under the coordinate chart corresponding to Boyer-Lindquist coordinates x a and their 
conjugate momenta pp by T^e,l z ,q for a given set of constants of motion. 

The basic topological properties of T^^e,l z ,q i n the case of a bound orbit of the first 
kind can be readily inferred from the equations determining the conjugate momenta, 

p t = -E, p^ = L Z) 

A 2 p 2 r = R, p 2 9 = Q. 1 ] 
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Firstly, since the orbital revolution progresses continuously in time, T^ e ,l z ,q is not 
compact in timelike directions. Secondly, the motion is bound both in the r- and in 
the ^-domain. Moreover, since p r is solely a function of r and, likewise, pe is a function 
of 9 only, both radial and polar motion are of compact support and of the libration 
type provided that the orbit is bound and stable. Otherwise, the orbit either plunges 
and, hence, the particle approaches the horizon towards future timelike infinity or, if the 
orbit is only marginally stable, it asymptotically approached the inner turning point. 
Thirdly, the azimuthal motion is a rotation and, of course, it is of compact support as 
well because = is to be identified with <fi = 2ir. 

According to a theorem by Arnold ||, that was formulated in the framework of 
non-relativistic mechanics, the phase space trajectory of a fully integrable system of N 
degrees of freedom is diffeomorphic to the iV-torus T N = (Sx) N whenever the motion of 
the system is of compact support. It seems plausible, that this theorem is applicable to 
relativistic problems as well if we account for the distinct role of the timelike coordinate 
which is related to the non-compactness of the worldline of a particle. Accordingly, as 
a manifestation of the weak equivalence principle, geodesic motion in Kerr spacetime 
has only three physical degrees of freedom because orbits are the same regardless of the 
mass (jl of the particle. Since T^ e ,l z ,q is defined with respect to a particular coordinate 
system, there is a well defined notion of projecting this set in the directions of the pair 
of conjugate variables (t,p t ). The resulting projected set, T^e,l z ,q, resides in the six- 
dimensional space spanned by the spacelike conjugate pairs (r,p r ), (9,p$) and {(j),p$) 
and corresponds to a phase space. Indeed, from the topological considerations in the 
previous paragraph, we see that T^ jE ,l z ,q is diffeomorphic to the three-torus T 3 . 

Utilising this toroidal topology, let us now define for each spatial component a closed 
curve Ck C T^^z^.q which circumscribes the fc-th torus and can be contracted into a 
single point with respect to the other tori. If ip is the coordinate chart on the cotangential 
bundle producing the Boyer-Lindquist coordinates and its conjugate momenta, then 
V> -1 (Cfc) is the pre-image of the k-th loop in T* (Ai^ eTT ) . Furthermore, if we note that 
the canonical one-form of Hamiltonian mechanics, © = —E dt + p r dr + pedO + L z d<p, 
is the one-form corresponding to the relativistic four-momentum p, i.e., 

(&,p)=PaP a = -fi 2 , (18) 
then we can define an action variable Jk for the k-th spatial component as 

J k = / 0. (19) 



2tt 



It might appear that such an action variable ultimately depends on the choice of the 
closed curve Ck which is defined with respect to a particular coordinate system and, 
thus, one would assume that it is not a geometric invariant. However, Jj~ is independent 
of the shape of Ck because is closed on 4>~ 1 (T^ e ,l z ,q), i-e., d0 = 0. The closedness of 
the canonical one-form is shown by substituting the expressions ( [TTD for the conjugate 
momenta and applying the rules of exterior differentiation^. The generalised Stokes 

§ It is actually the separability of the components of motion which renders the canonical one-form 
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theorem then yields the desired property of 0. Furthermore, since transformations 
between coordinate charts preserve topological features, we could have defined the curves 
Ck with respect to any coordinate system and then found their pre- images in T^A^Kerr)- 
It is only the topological features of Ck with respect to the three-torus that matters! 
Therefore we conclude that the action variables ([19]) are, indeed, geometric invariants. 
For convenience, however, we shall denote these variables as J r , Jg and as if they 
were specific for the Boyer-Lindquist coordinate system. 

Since motion in Kerr spacetime is separable in the coordinates r, 9 and 0, we may 
choose each Ck such that it is located completely within the (x fc ,pfc)-plane. Then the 
action variables can be easily calculated from cyclic integrals over the spatial conjugate 
momenta in the Boyer-Lindquist coordinate representation: 



As can be seen from the above expressions, J r , Jg and are functions of the mass jj 
and the three constants of motion E, L z and Q. 

The generator of the transformation to the action variables Jk is given by W = J ©. 
It is not globally defined, but changes by an amount of 2k Jk when integrating over the 
path Ck- As a consequence, if the spatial coordinate x k goes through N complete 
cycles and the other coordinates are kept unchanged, then the generalised coordinate 
w k = dW/dJk associated with J k changes by 27rA".[| This result implies that Wk is a 
phase or angle variable and, since Aw k = UkAT, the constant u>k can be interpreted as 
a frequency. Moreover, as a consequence of the geometrically invariant definition of 
the action variables Jk, the frequencies u>k specify fundamental properties of the orbital 
motion independent of any particular coordinate representation. Therefore, they are 
called the fundamental frequencies of the system. However, these frequencies are subject 
to a straightforward physical interpretation only with respect to the spatial coordinates 
r, 9 and 0, when the paths Ck are correspondingly chosen as cycles of, respectively, 
radial, polar and azimuthal motion. Accordingly, we shall use the notation u r , ug and 
ujrj) for the fundamental frequencies of an orbit in Kerr spacetime. 

The standard procedure of determining fundamental frequencies is to find the 
explicit form of the Hamiltonian in the action-angle representation, H^ & \ and to 
calculate the frequencies from the partial derivatives with respect to the action variables 




(20) 



(21) 



(22) 



dH (aa) 



(23) 



^ k = Wk 



being closed on the dynamical submanifold in the cotangential bundle, 
j Such a change has to be understood as a virtual displacement. 
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Unfortunately, the definition of the radial and polar action variables in terms of the 
non-trivial integrals fl20| ) and fl2~T|) does not admit an explicit inversion because none 
of these integrals can be solved analytically. However, as is shown in Appendix A , 
the derivatives / dJk and, hence, the frequencies can be found even without 

knowing the functional form of the Hamiltonian if( aa ) if the theorem on implicit 
functions is employed. This procedure results in equations ( |A.11| )-( |A.13 ) of Appendix 
[A| and the fundamental frequencies are therefore given by 



2tt 



2a 2 zl 



E{k) 



K{k) 

2it 4a 2 z + 4 X(n,r 2 



X(r 1 ,r 2 ) + 2Y(r 1 ,r 2 ), (24) 



^ p ^- E ^ + W + y^T) K{kl (25) 

2tt a 2 z\\K{k) - EffM + K(k)Y(r u r 2 ) 



fjw* L z [H(zi.,k) - K(k)]X(r h r 2 ) + K(k)Z(n,r 2 ) 

where r\ and r 2 are the turning points of radial motion, z± are the two roots of the 
equation Q(z) = when substituting cos^ = z in 0, k = z 2 _jz\, [3 2 = a 2 (fi 2 — E 2 ), 
K(k), E{k) and H(z^_,k) are, respectively, the complete elliptical integrals of the first, 
second and third kind, 

r /2 di> 

K(k) = / J. (27) 

</o VI - ft sin yj 

tt/2 



E(k) = / \Jl — k sin ip dip, (28) 
Jo 

U(zi, k) = f' 2 ^ } (29) 

Jo (l — z 2 _ sin 2 ip) \/l — k sin 2 ijj 

and X(ri,r 2 ), Y(ri,r 2 ) and Z(r\,r 2 ) are radial integrals defined by 

T 2 dr 

x(r " r ^L vn- (30) 

H 2 r 2 dr 

^=/ n W (31) 

f r2 r \L z r - 2M(L Z - aE)\ , . 

The above equations for the fundamental frequencies, however, have two 
deficiencies: Firstly, they are suggestive of being dependent on the mass /x of the 
particle which is, of course, not the case. For this reason, we will cast these equations into 
a dimensionless scale-invariant form, in which // does not appear anymore. Secondly, the 
radial functions X, Y and Z are not proper integrals because the integrated functions 
are divergent at the turning points r\ and r 2 . Fortunately, they are transformed into 
well-behaved integrals by virtue of the substitution r = Mr = pM/(l + ecosx), where 
X varies from to 2n as r goes through a complete cycle. Although this substitution 
originates from the pure Kepler ellipse, it is actually still expedient even for relativistic 
orbits that may be very different in shape. In analogy to non-relativistic orbits, p 
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is called the semi-latus rectum and e is the eccentricity of the orbit. Working out 
this substitution, yields the following equations for the dimensionless fundamental 
frequencies Uk = Mu>k' 

Ur ~ (l-eS)A' (33) 

»9 = (34) 

. (Z-L z X)K(k)+L z XU(z 2 _,k) . > 

^ = j , (35) 

where the tilde on top of each symbol indicates the corresponding dimensionless quantity 
(see Appendix ij| ) and the constant A is defined by 

A= (Y + a 2 z 2 + X)K(k) -a 2 z 2 + XE(k). (36) 

The dimensionless radial integrals in the ^-representation are given by 

i^^M^^ (37, 

^•T^-^r (1+ xw <38) 

Z(ri,r 2 )= / — ~ r— - , (39) 



!- e2 Jo H{x)^fW) 

where the functions G, H and J are defined by 

J( X ) = (1 - £ 2 )(1 - e 2 ) + 2 (l - £ 2 - l^pj (1 + ecosx) 

a 2 (i-£ 2 ) + L 2 + Q 



- 9 . 3 + 4 

+ <!(l-£ 2 )- --- + 

1 — e z p 

\2 



^ 2 



x (1 + ecosx) , (40) 

2 a 2 

H( X ) = 1 - -(1 + ecosx) + —(1 + ecosx) 2 , (41) 
p p z 

G( X ) = L, - 2 ^~~ aE \ i + ecosx). (42) 
p 

In the above equations, the constants of motion E, L z and Q are understood as functions 
of the orbital parameters p, e and z_ = cos 6L . Although analytical expressions for these 
functions can be derived, it is hardly feasible to substitute these expressions for E, L z 
and Q. For given orbital parameters, it is much easier to calculate the constants of 
motion numerically, as explained in [Appendix B| , and to substitute the resulting values 



into G(x), H{x) and J(x). 

Apart from the three frequencies u r , ujq and u^, a fourth constant of motion 
is obtained, which is associated with the timelike generalised coordinate X^Jt) = 
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XP^O) + 7T in the action-angle variable representation. According to equation ( A.lOj) 
of |Appendix A| , this constant is given by 

= [W{r u r 2 ) + fl 2 4£I(n,r 2 )] Kjk) - a 2 z 2 + EX( ri , r 2 )E(k) 
M7 [F(r:,f 2 ) + ^41(1^2)] if(fc) - ^^(r^r^lfc) ' ( ] 
where the function W(ri,r 2 ) is defined by 

W(r u r 2 ) = I" rlr* + a*Er-2a(L I -aE)}dr ^ 
J ri Av^R 
In dimensionless form, 7 can be expressed as 



7= A 



(W + ~a 2 z%EX)K(k) - a 2 z 2 + EXE(k) , (45) 



where 



with the generalised coordinate X? s is Pq = p t = —E, i.e., the same as in the Boyer- 
Lindquist representation. For this reason, one would expect that X,° > is in some sense 



W = — -^—W( )= r p2F{x) dx (46) 

M l-e 2 lri ' r2j J (1 + eMSxWVJW' 

and the function F is defined by 

F( X ) = E + ^(1 + ecosx) 2 - + ecos X ) 3 . (47) 

As for the physical interpretation of 7, we note that the conjugate momentum associated 

-0 

'(aa) 

mnn Ann ti r/-\i l I sA nxrr\ j^/tt r n nr V 

k (aa) 

an interval of coordinate time. However, it cannot be equal to t — to as equation (|14|) 
tells us that the relation between the elapse of proper time r and the corresponding 
interval of coordinate time, t — to, is not linear. On the other hand, there are three 
particular scales of proper time, r r , tq and 7^, associated with the three fundamental 
frequencies, u r , ujq and u^. Substituting each of these constants into the expression for 
X?^, three scales of coordinate time are produced: 

T k = X ( ° aa) (r fe ) - X ( ° aa) (0) = (48) 

This seems to be a natural definition of the coordinate-time scales associated with the 
fundamental frequencies. As we shall see in following Section, the above definition, 
indeed, reproduces the known periods of coordinate time in the limits of, respectively, 
circular and equatorial motion. Hence, we are led to the conjecture that 7 can be 
interpreted as a gravitational Lorentz factor which relates the three dynamical scales of 
proper time, r r , r e and r^, to the corresponding intervals of coordinate time, T r , T e and 

In the general case, this conjecture has the following implication: The ratio of time 
scales associated with two particular components of motion, say, radial and polar motion, 
is the same regardless of the frame of reference in which these time scales are defined. 
For example, if an observer at infinity saw the intervals of time T r and Tq, he would 
find that the ratio T r /Tg is equal to r r /rg. Of course, this is exactly what one would 
expect, provided that the fundamental frequencies truly arise from invariant properties 
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of orbital motion. In conclusion, the ratio of time scales should be independent of the 
frame of reference in use. With respect to coordinate time t, this is inherently ensured 
by the above definition of T r , Tq and T^. 

Ultimately, the significance of the fundamental frequencies lies in the possibility 
of expanding a dynamical quantity q{r) in a Fourier series with harmonics u>kim = 
ku r + Iujq + mu^ of these frequencies Q: 

q( T ) = ^2 qklm ex P( iuJ ki m r) = exp(ik ■ w), (49) 

k,l,m £ 

where k = (k, I, m) and k -w = kw r + lw e + mw^ = LOumJ according to definition of the 
action variables w r , w e and The Fourier coefficients q^ of the above expansion are 
given by 

j rlit r2iT />27T 

q^ = / dw r / dw e / du>^ q(k ■ w/u£) exp(— ik ■ w). (50) 

(2tt) Jo Jo Jo 

Therefore, dynamical properties of the particle are completely specified by discrete sets 

of numbers such as {q%}% & i?,- This is the very essence of bound geodesic orbital motion 

in a Kerr spacetime being quasi-periodic rather than quasi-chaotic, as was proposed due 

to the ergodic properties of these orbits. 



4. Limiting cases and numerical results 

As an important test to support the validity of the fundamental frequencies derived 
in the previous Section, we will consider two special cases of bound motion in Kerr 
spacetime: Firstly, the limit of vanishing eccentricity and, secondly, purely equatorial 
orbits. In particular, we will see that previously derived equations for the periods and 



the corresponding azimuthal angles of circular orbits |L0| can be reproduced from the 
general formulas. Briefly summarised, the features of orbits in the two limiting cases 
are as follows: 

• Circular orbits (e = 0): In Schwarzschild spacetime, Kepler's third law ([[]) holds 
exactly for the period of azimuthal revolution with respect to coordinate time, 

~3 /2 

i.e., = 2irMr . For circular orbits in Kerr spacetime, there are two distinct 
frequencies depending on whether the particle co-revolves or counter-revolves with 
the rotation of the black hole. If the motion is not confined to the equatorial 
plane, then the particle changes periodically its altitude such that 2n/uj0 is twice 
the elapse of proper time between two subsequent passages of the particle through 
the equatorial plane. The frequencies ujq and are generally incommensurate. 
As a consequence, the azimuthal angle $e = 2ttuj^/uj0 that is accumulated as the 
particle completes a cycle of polar motion is different from 2n. 

• Equatorial orbits (Q — 0): As a consequence of the different shape of the potential 
of radial motion as compared to the Newtonian potential, non-circular orbits do 
not close on themselves as the particle revolves around the gravitational centre. 
The total azimuthal angle $ r swept out by the particle as it moves from, say, the 
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apastron at radius r a — p/(l — e) to the periastron at radius r p = p/(l + e) and 
back to the apastron is given by $ r = 27ru) ( f > /u r . 

Actually, there is another limit, called polar orbits, for which L z = 0. In this case, the 
particle crosses the spin axis of the black hole and the inclination parameter is Z- — 1. 
However, it seems unlikely that polar orbits are astrophysically relevant and we only 
refer to [[□]] for a general discussion. 



4-1. Circular orbits 

Orbits of zero eccentricity, e = 0, are called circular. In this case, the particle moves on 
a spherical surface of constant radial coordinate, r(r) = pM = r , within the bounds 
of the inclination 9 given by < 6 < ix — Since the perihelion and aphelion radii 
coincide, we have r p = r a = ro and Rir§) = 0. In consequence, the integrals given 
by equations (PU|)-(P2D and ( pH| ) are undefined in the limit of vanishing eccentricity. 
The corresponding quantities in the ^-representation, however, have well-defined limits 
because the functions F, G, H and J are constant if e = 0. Hence, we define 

A-(r ,r ) = ^_ WM 



1/2' "V'U,<U, ^ jl/2 , 

(51) 



A ir 0^circ -"circ^circ 



circ 



,1/2' -V'U,'uy i/ 2 , 

A 4 '-'circ '"O-ncirc^circ 

where the constants F circ , G C i rc , #drc and J C j rc are obtained by substituting e = in 
equations (|40|)-(|42l) and fl47). The frequencies of the angular components of motion are 
therefore given by 

~ (circ) 

r ~~ (f 2 + ~a 2 z\ )K(k) - a 2 z 2 + E(k) ' K ' 

~(circ) _ 7Vf3z + 

We " 2[(r 2 + ~a 2 z\)K{k) - ~a 2 z\E{k)\ ' 1 ' 

~ (circ) ggjro - aL z )Ao ^(fc) + LJI(£, fc) 

^ (fg + 5 2 4)iT(A;) - a 2 z\E{k) ' 1 J 

The existence of a distinct radial frequency even in the case of circular motion 
indicates that the limit of the integrated azimuthal angle over a cycle of radial motion, 
$ r , does not become 2tt as the eccentricity approaches zero. Of course, this is well 
known from equatorial orbits. On the other hand, we can calculate the integrated 
azimuthal angle $0 as the particle goes through one cycle of polar motion by combining 
the equations of motion (|IlD and (0). The result is 

a(2Er — aL z 



$6 



(3z + 



A 



-K(k) + L z n(zt,k) 



(55) 



where Aq = r 2 . — 2ro + a 2 . When comparing this expression for $0 with the quotient of 
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cI>0 Circ ' ) and oj0 Cnc \ we see that the following well known equation holds: 

~ (circ) 

*' = 2 *%^y ( 56 ) 

Furthermore, calculating the period of coordinate time Tq for a cycle of polar motion 
from equation (|H|), yields 

T e = J- [M±3!lz^l^K(k) + ~a 2 z\E[K{k) ~ E(k)}\ . (57) 
(3z+ { Ao 1 J 

If the constant 7 given by equation fl45[) is also calculated in the limit e = 0, we find 
that 

T. = ^, (58) 

(circ) ' v ' 

which suggests that 7( circ ) can be interpreted as the Lorentz factor of circular motion in 
Kerr spacetime. 

4-2. Equatorial orbits 

If Q — 0, then the equation for the turning points of polar motion, Q(z) = 0, admits 
the solutions z_ = and z+ = 1 + L 2 z /f3 2 > 1. As a consequence, the motion of the 
particle is confined to the equatorial plane, i.e., 9(r) = tt/2. Since k = (z^/z + ) 2 = 0, 
the elliptical integrals become K(0) = E(0) = 11(0,0) = tt/2. The resulting expressions 
for the fundamental frequencies of equatorial orbits are thus given by 

^ = TT^V' (59) 

(1 - e 2 )Y eq 

" eq 

J«> = f*. (61) 

* eq 

As in the case of circular orbits, we still have three different frequencies. 

Combining the equations of motion ( |I0|) and ([12]), we can calculate the integrated 
azimuthal angle $ r between two successive passages through, say, the periastron and 
we find that 

2(1 - e 2 ) 

$ r = 2Z cq = -± J -Z cq . (62) 

P 

According to equation (|HD, the corresponding period of coordinate time is given by 
T r = 2W eq . Using 7^) = W eq /Y eq , we obtain 

T 2 ^ (6q) f R *\ 
" = ~Jf^~ ( } 

and one can see that 7^ acts as a Lorentz factor which relates the elapse of proper 
time over one cycle of equatorial motion with the corresponding period of coordinate 
time. 
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4-3. Keplerian orbits 

In the Newtonian limit, the constants of motion asymptotically approach the values of 
Keplerian orbits which are given by equations (|B.26| ) in |Appendix B| . Since the Keplerian 
values of energy and angular momentum, naturally, are independent of the spin of the 
black hole, the problem reduces to that one of motion in Schwarzschild spacetime and 
we can set z_ = 0. Thus, we obtain the the asymptotic expressions 

Fix) ^ 1 + OGtT 1 ), G( X ) * p 1 ' 2 + 0(p- 1/2 ), 

H( X ) z> 1 + Ofer 1 ), J( X ) * (1 - e 2 ) V 1 + 0(p' 2 ) 

which lead to 

W ~ Y, X^ Y ~ \ n , Z p^X, (65) 

1 — e z (1 — e z )U K 

where f^K is the Keplerian frequency given by equation (|Tj). 

Substitution of the above asymptotic expressions into equations (|33])-(|3"5"f) and 
results in 

7 ~ 1, Q r ~Q ~Qf~ il K . (66) 

The three components of motion therefore become commensurate in the Newtonian limit 
and the corresponding frequencies asymptotically degenerate. In the nearly Newtonian 
regime, where the three frequencies are close to Qk yet slightly different from each other, 
we expect that there will be a minute perihelion drift (due to the small difference between 
uj^ and uj r ) as well as a slow precession of the orbital plane (due to the difference between 
Utfy and ojq). In order to calculate this precession, an approximation of the fundamental 
frequencies including terms of next lower order in p would be required. However, this is 
left for future work. 



4-4- Numerical case study of a rapidly rotating black hole 

Let us now consider bound orbits around a nearly extreme Kerr black hole of spin 
a = O.998.0| Such a black hole is particularly interesting because effects due to spacetime 
dragging caused by the high angular momentum of the black hole are strong and stable 
orbits exist even very close to the horizon (which is located at r = rn ~ 1.0632M). 
Subsequently, we will explore the parameter space in all three dimensions and analyse 
the behaviour of the fundamental frequencies as each of the orbital parameters is varied. 
The numerical evaluation of equations d33|)-(p5|) and a general algorithm to calculate 
the constants of motion E, L z and Q for arbitrary orbital parameters p, e and 6- was 
implemented in MATHEMATICAL 

The chosen value of the spin is presumably a critical value as black holes which are accreting material 
tend to be buffered against further increase of spin if a w 0.998 

+ The procedures were tested, for instance, by comparison of large numerical samples to the results 
of Wilkins for circular orbits around an extreme Kerr black hole and those of Cutler et al jl4| for 
elliptical orbits in Schwarzschild spacetime. They are available from the author on request. 
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To begin with, we choose e = 1/3 (in this case the aphelion radius is twice the 
perihelion radius) and 6L = 7r/3 (corresponding to 30° of maximum altitude from the 
equatorial plane). The fundamental frequencies of the prograde and retrograde orbits 
as functions of p are shown in figure [I]. Basically, one can see that the frequencies u r , 
uj^ and ujq increasingly deviate from each other as the distance from the black hole 




1.6819 3 6 9.41628 30 100 

P 

Figure 1. Dimensionless fundamental frequencies u) r , ujq and as functions of p for 
orbits of eccentricity e = 1/3 and minimal inclination = 7r/3 to the spin axis of the 
black hole (long-dashed lines: cu r , short-dashed lines: u>g, solid lines: The set of 

lines diverging towards p » 9.416 corresponds to retrograde orbits (lower binding 
energy, smaller angular momentum), whereas the other set of lines which diverge 
towards p « 1.682 are the frequencies of prograde orbits (higher binding energy, larger 
angular momentum). The two vertical grid lines mark the locations of the marginally 
stable prograde and retrograde orbits. The widely spaced dashed straight line is the 
"asymptotic Keplerian branch" which is given by the power law J7k = (16-\/2/27)p -3 / 2 . 




Figure 2. Separatrix of a Kerr black hole of spin a = 0.998. The two surfaces specify 
the smallest radial distance r which can be reached by particles moving in a stable 
orbit of given eccentricity e and inclination parameter z_ without plunging into the 
black hole (a high-quality version of this graph can be obtained from the author). 



becomes smaller. For increasing p, on the other hand, the fundamental frequencies 
tend to degenerate and, in fact, they converge towards the Keplerian frequency Qk = 
(16a/2/27)p -3 / 2 for e = 1/3. Thus, we call the Keplerian power law for given eccentricity 
the asymptotic Keplerian branch of the fundamental frequencies. 

The values of p, for which the frequencies of, respectively, prograde and retrograde 
orbits diverge, are located exactly where the respective orbits become marginally stable. 
In this case, stable prograde orbits only exist if p > pml ~ 1.682 and stable retrograde 

(r) 

orbits only if p > p m i « 9.416. Hence, the corresponding minimal perihelion radii of 
stable prograde and retrograde orbits with e = 1/3 and 6L = 7r/3 are, respectively, 
r{, p) w 1.261M and rj, r) » 7.062M. Note that u r and u) -»• oo as p -> p ms . This is a 
consequence of the orbits becoming "zoom whirls" which means that the particle moves 
from the apastron inwards and then revolves many times along an almost circular spiral 
before it actually reaches the periastron.Q In the limit of a marginally stable orbit, the 

* Equatorial zoom whirl orbits have been recently investigated [|l5| . 
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periastron is approached asymptotically. The elapse of proper time as the particle moves 
from apastron to periastron, thus, ever increases as p — > p ms and the radial frequency 
approaches zero. If the particle represents an astrophysical object, however, radiation 
reaction will finally act on a time scale that is comparable to the dynamical time scale of 
radial motion (~ 2n^/uj r ) and, as a consequence, it will drive the particle within finite 
time over the stability threshold so that it plunges into the black hole. 

In parameter space, the locus of all points corresponding to marginally stable orbits, 
P = Pms(e, is called the separatrix. A plot of the separatrix for a = 0.998 is shown 



0_ = jt/3, p = 12 6L =tt/180, p = 12 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



e e 



Figure 3. Dimensionless fundamental frequencies u) r , Cjg and cj^ as functions of the 
eccentricity e for different values of p and 6L (long-dashed lines: Q r , short-dashed 
lines: Qg, solid lines: ui^; the frequencies of retrograde orbits enclose those of prograde 
orbits). In the bottom panels, only the frequencies of prograde orbits are plotted since 
there are no stable retrograde orbits if p = 6. 
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in figure 0, where the perihelion radius r p = p/(l + e) is used as parameter in place of 
p. There are two sheets, one corresponding to prograde marginally stable orbits and 
the other one corresponding to retrograde marginally stable orbits. The two sheets 
smoothly join at z = 1, i.e., where the polar orbits with L z = are located. As 
Z- = cos 6L becomes smaller, however, there is an increasing space between the two 
sheets, in which only prograde but no retrograde orbits reside. As the eccentricity 
increases, both prograde and retrograde orbits exist closer to the horizon due to the 



e = 1/3, p = 12 



e = 9/11,p = 12 




0.03 

0.02 
0.015 

0.01 




0.2 0.4 0.6 Oi 

Z- 




Figure 4. Dimcnsionlcss fundamental frequencies Co r , &g and uj^ as functions of the 
minimal inclination to the spin axis of the black hole, for different values of e and 
p. (long-dashed lines: H) r , short-dashed lines: Qg, solid lines: cD^; the frequencies of 
retrograde orbits enclose those of prograde orbits). In the bottom panels, only the 
frequencies of prograde orbits are plotted since stable retrograde orbits with p = 4 
do not exist at all. Even the prograde orbits are only stable if the altitude from the 
equatorial plane is lower than roughly 37r/8 for e = 1/3 and about ir/3 for e = 9/11. 
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higher momentum of the particle at the periastron as compared to circular orbits of the 
same radius. In particular, there are prograde orbits of eccentricity e > 0.5 with z < 0.5 
for which the periastron is very close to the horizon (within 5r ~ 0.1). 
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Figure 5. Scales of coordinate time T r , Tg and in hours as functions of the semi- 
major axis a = pM / (1 — e 2 ) in astronomical units (1 AU = 1.496 • 10 13 cm) for a central 
black hole of mass M = 3- 10 6 Af so i (long-dashed lines: T r , short-dashed lines: Tg, solid 
lines: T^). As in figure |l|, the eccentricity of the orbits is e = 1/3 and the inclination 
parameter 6- = n/3. The set of lines diverging towards a w 0.3137 AU corresponds to 
retrograde orbits, whereas the other set of lines which diverge towards a « 0.0560 AU 
are the frequencies of prograde orbits. The two vertical grid lines mark the locations 
of the marginally stable prograde and retrograde orbits. The widely spaced dashed 
straight line is the Keplerian orbital period given by = 4ir 2 a 3 / 'GM (in physical 
units). 

Figure |3] shows the change of the fundamental frequencies with eccentricity for 
various choices of 6L and p. The main feature of these plots is the convergence of 
the frequencies to a single value as e — > 1. This simply reflects the fact that orbits 
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of eccentricity very close to unity are almost unbound and the apastron is located far 
away from the region of strong gravity. In consequence, the particle is mostly moving 
in the nearly Newtonian regions of spacetime and the fundamental frequencies tend to 
degenerate. In the plot showing the frequencies of nearly polar orbits with Q_ = 7r/180 
and p — 6, the frequencies terminate just below e = 0.5 as prograde orbits of higher 
eccentricity do not exist if z ~ 1 (note that the value of the perihelion radius depends on 
e according to f p = 6/(1 + e) when comparing to the plot of the separatrix in figure §). 

The dependence of the fundamental frequencies on the inclination parameter 6L 
(which is the smallest angle of inclination to the spin axis of the black hole that can be 
reached by the particle) is shown in figure f| for different values of e and p. Apparently, 
both the radial and polar frequencies of prograde and retrograde orbits converge towards 
each other as 9- — > 0, i.e., in the polar limit. On the other hand, there is a crossing-over 
of the azimuthal frequencies at some finite value of 9- . As a consequence, the azimuthal 
motion of nearly polar retrograde orbits is faster than that of the corresponding prograde 
orbits and, thus, it seems that nearly polar orbits are markedly different from nearly 
equatorial orbits. We have to keep in mind, however, that the revolution of the particle 
around the black hole in the case 9- ~ is rather specified by the polar frequency ug, 
whereas determines an azimuthal precession around the spin axis of the black hole. 

Finally, in figure ^ the scales of coordinate time, T r , T e and T^, defined by 
equation ^ are shown in physical units as functions of the semi-major axis a for orbits 
of eccentricity e = 1/3 and inclination parameter 6L = ir/3. The mass of the black 
hole is chosen to be M = 3 • 10 6 M so i which is about the mass determined for the black 
hole in the centre of our Galaxy [jnj. It is interesting to note that for virtually all 
prograde orbits the time scale of polar motion T$ is very close to the the Keplerian 
orbital period Pk = 27r/f2K- However, this is only accidentally the case for the chosen 
set of parameters, as for lower eccentricity is above Pk, whereas it tends to become 
significantly smaller than Pk at higher eccentricity. 

5. Conclusion 

The investigation of bound geodesic orbits in Kerr spacetime presented in this article 
clearly illustrates that the properties of these orbits in the regime of strong gravity 
are profoundly different from Keplerian orbits in the Newtonian regime. In particular, 
there are three different frequencies of orbital motion, u r , ujq and u^, which increasingly 
deviate from the Keplerian frequency as the size of the orbit decreases and finally 
approaches the limit of a marginally stable orbit. This is a consequence of the radial, 
polar and azimuthal components of motion being incommensurate in Kerr spacetime. 
As a consequence, the shapes of orbits in the vicinity of a Kerr black hole are exceedingly 
complex. 

The explicit formulas for the fundamental frequencies u r , ujq and ui^ presented in 
Section § are exact in the test particle limit and useful approximations if radiation 
reaction is adiabatic, i.e., if the time scales associated with the frequencies are small 
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compared to the time scale of radiation reaction. Based on the constant of motion 7 
which is associated with a timelike generalised coordinate, we have proposed a conversion 
of the invariant fundamental frequencies into coordinate-time scales, T r , and Tg, 
which should be relevant to an observer at spatial infinity and, thus, gravitational wave 
astronomy. Once these time scales are established, one can compute instantaneous 
gravitational wave forms in the adiabatic approximation. In this context, the term 
"instantaneous" means an interval of time, over which no significant changes of the 
constants of motion are accumulated due to radiation reaction. From the resulting 
spectrum, an estimation of the power emitted in the frequency band of a detector could 
be made, without making use of an advanced theory of radiation reaction. 
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Appendix A. Determination of the fundamental frequencies 



Let P, 



(aa) 



p(aa) _ w h ere j k are the action variables defined by equations 



/ (aa) ( 



—jj, /2,E,L z ,Q) be the momenta given by P 



(aa) 



Pt 



■ —E and 

/, - •.//,.. wijlvjv 1 jl j v -./ aiv ilia av niiii vaamiiiv urumu in vi[uai juiin \ _ v i / \ ]22| ) . If We 

denote the Jacobian matrix of / by Df, then, by the theorem on implicit functions, 
Df ■ D(f~ 1 ) = Df ■ (Df) -1 = I, provided that / is non-zero and the Jacobian 
does not vanish \\L7\ . Since —fi 2 /2 is the invariant value of the Hamiltonian, we can 
substitute —fi 2 /2 = ij( aa )(— E, J^), where is the Hamiltonian in the action-angle 

variable representation. For brevity, let us subsequently use the symbol H to denote 
the Hamiltonian in that representation. Moreover, two rows of the Jacobian matrix are 



trivial due to the identities Pq^ = —E and = L z . Thus, the equation Df-D(f~ 



I 



reads 



/0 
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dH 
dJe 
dH 



-1 
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dE 
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8L Z 
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8Q 
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8Q 
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dQ 


dQ 


dQ 


dQ 


dE 
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dJg 


dJj, 



(a.i; 



The above matrix equation can be split into four non-trivial sets of linear equations in 
the eight unknowns — 



~fi and H : 




2 r r 2 dr 

J r 



>ri VR 



(r 2 +a 2 )P 



+ a(L z — aE) 



2a*Ef t 



t/ 2 cos 2 6>d6> 
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(A.4) 
(A.5) 



r, r f t/ 2 cot 2 8d8 

where P = (r 2 + a 2 )-E — aL z , and the coefficient matrix A is given by 

rr 2 r 2 dr f"2 dr 



9 f r 2 r^dr _ rrn dr_ 

Jn Jri Vr 

o„2 f 71 "/ 2 cos 2 8 d9 d8 



_ , (A.6) 

Using the definitions and ( fP| ) for the radial integrals and the identities 

-/2 de x 



~ f3z- 



-K(k), (A.7) 



w/2 cos 2 9 

C -^-d6= Z -±[K(k)-E(k)], (A.8) 

n/2 cot 2 (9 z^ 

C -^-d6= Z -±[U(z 2 _,k)-K(k)}, (A.9) 

where the elliptical integrals K(k), E{k) and IT(— z 2 _ , k) are defined by equations (f27j)- 
(p9j), the solutions of the above systems of equations for — and ^j- are given by 



dH _ K(k)W(r u r 2 ) + a 2 z\E [K(k) - E(k)] X(r u r 2 ) 
~dE ~ K(k)Y(n,r 2 ) + a 2 z\ \K{k) - E(k)\ X(n,r 2 ) ' 
dH _ nK(k) 
d~T r ~ K(k)Y(n, r 2 ) + a 2 z 2 + [K(k) - E{k)} X(n, r 2 ) ' 

dH_ nPz + X( ri ,r 2 ) 

dJ e ~ 2{K(k)Y(r u r 2 ) + a 2 z\ [K(k) - E{k)} X(r u r 2 )} ' 
dH K(k)Z(r u r 2 ) + L z [U(z 2 _,k) - K(k)]X(r u r 2 ) 



(A.10) 
(A.ll) 
(A.12) 
(A.13) 



dJ^ K(k)Y{r u r 2 ) + a 2 z\[K(k) - E(k)]X(r 1 , r 2 ) ' 

The validity of the above solution depends on the Jacobian matrix Df being non- 
singular. In fact, det Df ^ if and only if 



1 _ E ( k ) 



dH 

dJ r 



2 



t^O, (A.14) 



K(k) 

as is easily shown by using the second equation of system ( |A.3| ). The above condition 
is fulfilled whenever k > 0, i.e., if the orbit is non-equatorial. Since Q = for 
equatorial orbits, this is immediately clear from the above systems of equations as 
well. Nevertheless, the frequencies of equatorial orbits can be calculated by taking the 
limit k — > of the general expressions for k > and, as demonstrated in Section |4.2j , 
one finds agreement with the results obtained from a direct calculation based on the 
equations of motion. 
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Appendix B. Calculation of the constants of motion 



In order to find the constants of motion E, L z and Q for given orbital parameters p, e 
and 6L, one has to solve the following set of equations: 

^ = <=► i?(r)=0, (B.l) 

6(0) = 0, (B.2) 



where i?(r) and 0(0) are given by equations (0) and (H). The roots of these equations 
correspond to the turning points of, respectively, radial and polar motion, i.e., we can 
substitute r = r p = pM/(l + e) or r = r a = pM/(l — e) in the first equation and 
9 = 6- in the second equation. Thus, a system of three equations is obtained which 
can be solved for the three unknown constants of motion, provided that e ^ 0. In the 
case e = 0, on the other hand, we have got to make use of the supplemental constraint 
R'(r ) = 0, in addition to R(r ) = and 6(6L) = 0, because circular orbits occur if 
both R(r) and its radial gradient R'(r) vanish at the same point, r = r Q . 
Introducing the dimensionless quantities 

L z — tti Q — — 



a ~_E 



(B.3) 



fiM' /x 2 M 2 ' 

we can actually set up a completely scale-invariant formulation of the problem. Hence, 
it is sufficient to calculate E, L z and Q only once for given orbital parameters. The 
corresponding physical values for a particular black hole of mass M and a particle of 
mass [x can then be inferred from the above definitions. 

To begin with, we re-arrange 0(0_) = in order to express Carter's constant Q in 
terms of 6_, E and L ? , 



Q 



a 2 (l 



E') + 



(B.4) 



where z_ = cos6L. If this expression is substituted for Q in the equation R(r) 
Li- 2 M~ 4 R(r) = 0, we obtain 



where 



R{f) = f{r)E 2 - 2g{r)EL z - h(r)L 2 z - d(r), 



f(r)=f i + a 2 f(f + 2) + z 2 _A 
g(r) = 2af, 



h(r) = r(r — 2) + 



-A, 



1-z 2 ' 
d(r) = (f 2 + fiV)A, 

and A = f 2 — 2f + a 2 . Substitution of Q in R'(r) = yields 
— = f(r)E 2 - 2g'(r)EL z - h'(r)L 2 z - d\r), 



(B.5) 

(B.6) 
(B.7) 

(B.8) 
(B.9) 

(B.10) 
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(B.H) 
(B.12) 

(B.13) 
(B.14) 



d'{r) =2(2r- 3) r 2 + 2a 2 [(1 + z 2 _)r - z 2 ] . 

Thus, we shall define the following coefficients in terms of the orbital parameters: 

(f(r P ),g(r P ), h(f p ), d(f p )) if e > 0, 
(f(fo),g(fo),h(f ),d(fo)) ife = 0, 



(f2,g2,h 2 ,d 2 



(f{r a ),g(r a ), h(f a ), d(f a )) if e > 0, 
(f(r ),g>(r ),h'(r ),d'(r )) if e = 0. 



(B.15) 
(B.16) 



Eliminating L z from the set of quadratic equations, 



fiE 2 - 2 9l EL z - h x Ll - d x 



0. 
0. 



f 2 E 2 - 2g 2 EL z - h 2 L 2 z - d 2 

a quadratic equation in E 2 is obtained, 

(p 2 + Aria)E 4 - 2(k P + 2ea)E 2 + k 2 = 0, 
where e, rj, K, p and a are defined by the following 2x2 determinants: 
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The roots of equation (|B.18|) are given by 

■ n2 np + 2ecr ± 2y 1 o(oe 2 + pen — nn 2 ) 



E. 



± 



p 2 + Ar/a 



(B.17) 
(B.18) 

(B.19) 
(B.20) 
(B.21) 

(B.22) 



Actually, these roots can be calculated regardless of whether the corresponding orbits 
are of the first or the second kind, i.e., if they are stable or plunging. There is a 
simple criterion, however, which can be used to determine the kind of the orbit: Let the 
potential of the radial component of motion be defined by 



V 2 . 



(B.23) 



According to equation (0), V = (E 2 — R/Y, 2 ) 1 / 2 and the radial domain of motion is 
subject to the constraint V < E. Hence, the orbit is stable if 

dV 



(B.24) 
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If the gradient of the potential vanishes at r = r p , the potential has a local maximum 
at the perihelion radius and, in that case, the orbit is marginally stable. 

Substituting ±y/E± into the system ( |B.17| ) and solving for L z , yields four 



roots which satisfy both equations. Therefore, we obtain four possible solutions for 
the constants of motion with given orbital parameters, namely, (— — Li , Q^), 
(-E( p \ -Li p) , Q(p)), (E( p \ Li p) , Q(p)) and L { *\ QW), where E® = E. is the lower 
and E^ = E + the higher energy given by equation (|B.22|) . Numerical evaluation shows 
that L z p) > L z t) . In the first case, the particle has higher binding energy and co-revolves 
with the rotation of the black hole. Such orbits are called prograde. In the second case, 
the particle has lower binding energy and usually counter-revolves, except for some orbits 
close to the horizon of rapidly rotating black holes. These orbits are called retrograde. 
The solutions with negative energies correspond to the respective motions under time 
reversal. 

Finally, we consider the Newtonian limits of the constants of motion which are 
asymptotically approached as p — > oo. Using the formulas for the potential energy and 
angular momentum of Keplerian orbits in geometric units ||, 

U = - 1 —^, L=pW (B.25) 
2p 

and noting that L z = L cos(7r/2 — 9 J) = L sin 9_ = L(l — z 2 ) 1 / 2 , we obtain the following 
asymptotic expressions for the constants of motion:|| 

1 - e 2 

1 - E 2 ~ , L 2 ~ (1 - z 2 _ )p, Q ~ z 2 _p. (B.26) 

P 

Since Q + L 2 ~ L 2 , Carter's constant in the Newtonian limit is the squared angular 
momentum component parallel to the equatorial plane. 
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